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Abstract 
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I. INTRODUCTION 



Flavor changing interactions in general, and of top quark in particular, are a promising 
test ground for new physics. In the Standard Model (SM) processes such as t — > cg,'j,Z 
are absent at tree-level and highly suppressed by the Glashow-Iliopoulos-Maiani (GIM) 
mechanism at one loop. The branching ratios predicted in the SM jl| are of order of 
10^^°; 10^^^; 10^^'^ — 10^^^ for the decays t cg^cj^cZ and thus far from present and 
future reaches of e^e~ colliders or the LHC Models beyond the SM predict branching ratios 
which are orders of magnitude larger than this, and thus an experimental signal can be 
interpreted as a signal for new physics. An added bonus of top quark phenomenology is that 
it provides a window into the electroweak breaking mechanism, since the top mass is close 
to the scale of symmetry breaking. 

Previous analyses have explored several possibilities of enhancing flavor violating decays, 
mostly within supersymmetry. t —>■ cV {V = g, 7, Z) have also been explored in two Higgs 
doublet models 1111^, technicolor models Js^, top-color assisted technicolor models Q, models 
with extra vector singlets [5| and supersymmetry with llll . [1^ , and without 

[13^ . R parity. While supersymmetry appears to provide the most convincing scenario for 
new physics, most analyses concentrate on the minimal supersymmetric standard model 
(MSSM). Analyses have been performed initially in the constrained MSSM, where flavor 
changing interactions are driven by flavor universal supersymmetry breaking soft terms 
including left, right and intergenerational scalar quark mixing. Li et al. evaluated one- loop 
SUSY QCD and EW contributions and obtained branching ratios of 10~^ for t —* eg 
and 10"*^ for t —>■ cy, Z. Couture et al. included the left-handed squark mixing in their 
calculation [?! and obtained branching ratios of 10"^ for the gluon, 10~^ for the photon 
and 10~^ for the Z boson. Inclusion of the right-handed squark mixing 8], did not change 
branching ratios significantly. Lopez et al. further refined the calculations in the flavor 
diagonal case by including contributions from the neutralino quark-squark loops; in this 
case the branching ratios obtained were 10~^ for t ^ eg and 10"'' for t — 07, Z. In the 
non-universal case 
for t — i> C7, Z. 

In the unconstrained MSSM, assumptions about universality of soft supersymmetry 
breaking terms and new sources of flavor violation are included in the scalar quark ma- 
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trices. Liu et al. looked at the fully unconstrained MSSM and found decay branching ratios 
which can reach 10~^, 10~^ for t —>■ eg, and t —>■ cj, Z respectively; while Delepine and Khalil 
[l^ found, for the unconstrained MSSM with a light stop, branching ratios of 10~^ for the 
gluon decay and 10^^ for the photon. Note that not all authors use the same values for 
relevant parameters. Experimental limits on SUSY masses have decome more stringent, in 
particular, the gluino mass is taken to be higher in later papers, which effects the estimated 
branching ratios. 

However, not much work has been done on flavor changing interactions in the top 
quark decays in scenarios beyond the MSSM. We present here such an analysis, based 
on extending the SM group SU{3)c x SU{2)l x f/(l)y to a left-right symmetric group 
5f/(3)c X SU{2)l X SU{2)r x U{1)b-l (LRSUSY). This gauge symmetry allows for the 
seesaw mechanism within a super symmetric scenario and predicts neutrino masses and mix- 
ing naturally. Flavor violation in h decays have shown possible enhancement of results over 
the the MSSM ones: in particular restrictions on intergenerational left-left (LL), right-right 
(RR), and left-right (LR) mixing of scalar quarks is more restrictive than in the case of 



MSSM \U\. 

In this work, we present a full investigation of the flavor changing two body decays of 
the top quark in the LRSUSY. We investigate both the constrained case, in which the only 
source of flavor violation comes from the Cabibbo-Kobayashi-Maskawa (CKM) matrix in the 
quark sector (with the proviso that we assume it to be the same for left- and right-handed 
quarks), as well as in the unconstrained model, in which soft symmetry breaking parameters 
are allowed to induce flavor-dependent mixing in the squark mass matrix. Such a flavor 
dependent scenario is motivated by neutrino oscillations, which indicate the presence of two 
large mixing angles (between the first and second, and second and third generations). If 
LRSUSY occurs in nature as an intermediate symmetry in breaking down from a SUSY 
GUT scenario, such as 5*0(10), large neutrino mixing angles in the lepton mass matrices 
may appear not in the down quark mass matrices (because the right-handed charged current 
interaction is broken down at the unification scale), but in the squark mass matrices jlfil ]. 
The new mixing might affect h [16] as well as t flavor physics. 

Previous studies of unconstrained supersymmetric models have been based on either 

n nn 

the mass insertion 171 or on the mass eigenstate method 16, la^. In the mass insertion 

n 

framework [12l| one chooses a basis for fermion and sfermion states in which all the couplings 
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of these particles to neutral gauginos are flavor diagonal. Flavor changes in the squark sector 
arise from the non-diagonality of the squark propagators. Off-diagonal elements mix squark 

ri 

flavors for both left- and right-handed squarks. In the mass eigenstate method [18], squark 
mass matrices are given in the super-CKM basis, and are diagonalized by rotating the 
superflelds. In this basis, the up-squark and down-squark mass matrices are correlated by 
this rotation and thus not independent. Potential new sources of flavor violation arise from 
couplings of quarks and squarks to gauginos. This method has the advantage that, when the 
off-diagonal elements in the squark mass matrices become large, the method is still valid, 
unlike the mass insertion which is a perturbation-based expansion. We will return to these 
considerations when performing specific evaluations, and chose the mass eigenstate method 
for its greater flexibility. 

Our paper is organized as follows: In Section ^1 we present the general framework of 
the left-right supersymmetric model. The squark sector is then discussed in Section IIIII in 
both the flavor diagonal and non-diagonal scenarios, concentrating especially on the flavor 
changing mixing between the second and third generations. In Section IIV[ we give the 
full analytical expressions contributing to t cV {V = g, 7, Z) decays by providing 
expressions for the gluino, chargino and neutralino contributions separately. The numerical 
analysis of the decays is given in Section IVl where we concentrate on individual and relative 
contributions of gluinos, charginos and neutralinos to the branching ratios of the decays. 
The discussion is carried out in the flavor diagonal and non-diagonal scenarios. Section IVTl 
is devoted to our summary and conclusion. The chargino and neutralino mass matrices and 
their diagonalization procedures are summarized in Appendix El Appendix El includes the 
relevant Feynman rules. 



II. DESCRIPTION OF THE LRSUSY MODEL 



The minimal supersymmetric left-right model is based on the gauge group SU{3)c x 
SU{2)l X SU{2)r X U{1)b-l 03- The matter flelds of this model consist of three 
families of quark and lepton chiral superflelds with the following transformations under the 
gauge group: 



(3,2,1,1/3), Q'r 



(3,1,2,-1/3), 



u 



R 



Vr 



LI = I M ~ (1,2,1,-1), Vr=\ \ ~ (1,1,2,1), (2.1) 



where the numbers in the brackets represent the quantum numbers under S'f/(3)c x SU{2) l x 
SU{2)£i X U{1)b-l- The Higgs sector consists of bidoublet and triplet Higgs superfields: 

I I -(1,2,2,0), {z = u,d), 

Al=| ^ U (1,3,1,-2), 5,= \^ V I ~(1'3'1'2), 



A«=| ^_ ~ (1,1,3,-2), Sn=[^^ ^^,1^(1,1,3,2). (2.2) 

The bidoublet Higgs superfields appear in all LRSUSY to implement the SU{2)l x U{1)y 



symmetry breaking and to generate a CKM mixing matrix. Supplementary Higgs represen- 
tations are needed to break left-right symmetry spontaneously: either doublets or triplets 
would achieve this, but triplet Higgs bosons and A/j are chosen to facilitate the seesaw 
mechanism |21|. In supersymmetry, additional triplet superfields 6r are introduced to 
cancel triangle gauge anomalies in the fermionic sector. The most general superpotential 
involving these superfields in LRSUSY is: 

W = Y^^^Ql^,tT2QR + Y^l^Ll^,tT2Ln + YLR{LltT26LLL + LltT2ARLn) 

+fiLR [Tr(Ai5L + Ar6r)] + ^iijTr{iT2^JiT2^j) + Wmr (2.3) 

where Yg and Y^ are the Yukawa couplings for the quarks and leptons, respectively and Y^^ 
is the coupling for the triplet Higgs bosons. The parameters /i^ and are the Higgs mass 
parameters. Left right symmetry requires all Y-matrices to be Hermitean in the generation 
space and Ylr matrix to be symmetric. Here Wnr denotes (possible) non-renormalizable 
terms arising from higher scale physics or Planck scale effects j2^. The presence of these 
terms insures that, when the SUSY breaking scale is above Mv^^,, the ground state is R- 
parity conserving. In addition, the Lagrangian also includes soft supersymmetry breaking 
terms as well as F- and D-terms: 



C 



soft 



P^qY^^QI^^it^Qr + A^Y fll^^ir^LR + AlrY LR^Liir^dLLL + L'^it^ArLr) 
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+ 



[m, 



+ H.c] - [5/.,,$,$, + H.c] + [{m'QjijQliQLj + {ml^)ijQRi9R3\ (2-4) 

where Aq, Al^Alr are trilinear scalar couplings. The LR symmetry is broken spontaneously 
to t/(l)em through non-zero vacuum expectation values (VEV's) of the Higgs fields. These 
values are: 




d 





Ha 

0/ U^, 




where we already set to zero the CP-violating phase in the mixing of Wl and Wr. The 
non-zero Higgs VEV's break both parity and SU{2)r. In the first stage of breaking, the 
right-handed gauge bosons, Wr and Zr acquire masses proportional to i'a^, '^'5^ and become 
much heavier than the SM (left-handed) gauge bosons Wl and Zl, which in turn pick up 
masses proportional to and Kd at the second stage of breaking. 

Fermionic partners of gauge and Higgs bosons mix. In LRSUSY there are six singly- 
charged charginos, corresponding to A^. A/j, (p^, (t)d, ^l, and A^. The model also has eleven 
neutrahnos, corresponding to Xz, \z', ^v, 4>iu^ 4>2w 4>id^ 4>2di A?,, A^j, 5^^, and 5%^. The 
doubly charged Higgs and Higgsinos do not affect quark phenomenology, but the neutral 
and singly charged components do, through mixings in the chargino and neutralino mass 
matrices. 

The supcrsymmctric sources of flavor violation of interest in t decays in the LRSUSY 
model come from either the Yukawa potential or the trilinear scalar coupling. 

The interaction of fermions with scalar (Higgs) fields relevant for the quark sector has 
the following form 

Cy = ^qQl^uQr + ^qQL^dQR + H.c. (2.5) 

We present below a detailed analysis of fiavor violation in the model in both the flavor- 
diagonal and non-diagonal case, before proceeding with calculation of the branching ratio 
of t cV. 
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III. FLAVOR CHANGING IN THE LRSUSY 



When a basis, called the super-CKM, in which the quark states are diagonal, is used to 
express the squark mass matrix, flavor non-diagonal entries naturally arise. So, our analysis 
is focused on the effects of such scenarios in rare top quark decays in comparison with the 
flavor diagonal case. 

In the interaction (flavor) basis, ((5i^,<5/j), the squared-mass matrix for squarks is 

^ ( '^}LL + ^fLL + DfLL i^) lr) + Ff LR \ 

The F-terms are diagonal in the flavor space, FfL^jjui = m^, {FdLR)ij = jJ^ij^di tan/3)ljj, 
{FuLR)ij = fJ'i^Ui cot where tan/? = i^u/i^d is defined. The D-terms are also flavor- 

diagonal 

DfLL = m|cos2/3(T3/ - Q/sin^6'w)l3x3, 

DfRR = m|cos2/3Q/sin^6'wl3x3- (3.2) 

In the universal case, one has {'fn^f llj RR)ij — "^1^^*^*^' ^}lr — ^*f^f- To reduce 
the number of free parameters, we also consider the following parameters to be univer- 
sal, = ^susy'^ji' ^d,ij = ^^ij and Au^ij = A5ij. Note that due to the in- 
variance under both SU{2)l and SU{2)r, '^MLL.uijij the up sector can not be speci- 
fied independently from the matrix frtdLLdRRi down sector. They are related as 



2 _ 1^2 \ i^] 

u LL,u RR 



KcKM {mlLL,dRR) -^CKM' where Kqkm is the CKM matrix.^ 



The mass(-squared) matrix in the universal case for the U-type squarks then reduces to, 
in block form 



M 



Msusv + ^l, + (^3« - Qu sin^ 9w) cos 2(3 ( A + /i cot /3) 



2 _ I susy Uk 
Uk ~ 



Vfluk + ^^ cot 13) ^Insy + + slu^ 6'w COS 2(3 

and, for the D-type squarks, to 

^2 _ I M^^^y + m^^ +m|(T3d - (5dsin^6'w)cos2/? m^JA + /itan/?) 



rud^, (v4 + /i tan /5) + ^2 _^ rn%Qd sin^ 6'w cos 2(3 



^ For simplicity, we assume ^^^ckm ~ -^ckm which is a conservative choice and doesn't require new mixing 
angles in the right-handed quark matrices |l4|. 
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The corresponding mass eigenstates are defined as 




(3.3) 



QR 



where Tq^ q/j are 6x3 mixing matrices and g is a 6 x 1 column vector. 

These up- and down squark mass matrices are 6x6, but are written above in 2 x 2 block 
form. In the flavor diagonal scenario (constrained LRSUSY), each of these blocks has no 
non-zero off-diagonal elements. That is, there is no intergenerational mixings for squarks 
and the only source of flavor mixing comes from the CKM matrix. 

The way to induce off-diagonal entries radiatively is to consider the evolution of the 
squark and quark masses from the SUSY-breaking scale (where both quark and squark 
mass matrices are flavor diagonal in the same basis), down to the electroweak scale via 
renormalization group equations. Experimental bounds from — and — data 
I23I 1 involving the first generation are tightly constrained. The mixings between the second 
and third generations are, on the other hand, free j^]. So, in our analysis we assume 
significant mixing between the second and third generations in the up- and down-squark 
mass matrices and neglect those involving the first generation. 

The up-squark mass squared matrix in the {ul, cl, t^, ur, cr, tji) basis can be written as 



/l 





M 



2,FC 



M: 



susy 





1 





(7^)23 











\ 



{d^r""),, 



U )22 W(/ )23 
urn 1 (5^«)32 (5 

10 

U )22 {S§^)23 1 

V (5^^)32 {S^%3 {S§%2 1 



(5^-^^ 



{6^"^^ 



LR\ 

u J33 



iS§%3 



(3.4) 



Here we assume 



where (5^^)22 = m,(A + /icot/3)/MLy, (5^^)33 = mt{A + fi cot (3)/ M^,^,^. 
that all diagonal elements on the main diagonal are set to a common value M^^^y ^ and the 
flavor changing off-diagonal entries of each block are furthermore scaled with M^^gy to make 



^ We comment on the effect of relaxing this condition. Basically, there are non-negligible corrections coming 
to the 3rd and 6th diagonal entries of both the up- and down-squark matrices. The term to be added 
to unity in the up-sector is m^/M'^^^y, and in the down-sector ml/M^^^y with a contribution from the 
D-term. So, only in the up sector, the Mgugy ~ rrit case can give sizable contributions. However, in the 
flavor non-diagonal scenario, we don't consider such small A/gusy values. 
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the parameters dimensionless. Thus, these off-diagonal parameters are defined as 



_ (^M Ll)u (xLR\ _ ("^M LR)ij 



susy susy 
^2 



(O.- = ^if^ ^ = ^^^f^ , ^^J = 2,3. (3.5) 



T2 ' ^"U )ij ^ 

susy susy 

Here we still take the diagonal elements of muLR,uRL same as in the flavor diagonal case. 
The down sector case is completely analogous to the up sector. 

As mentioned in the Introduction, there are two approaches to compute the effects of 
the flavor changing parameters 5's to physical quantities. The first historical one is the 
mass insertion formalis m 117 1 and the second and more recent one is the general mass 
„e fonnal.. U^l. the .a. in.e.,o. fo^aUs.. the . te™s .present .ixing 
between chirality states of different squarks, and it is possible to compute the contributions 
of the first order flavor changing mass insertions perturbatively if one assumes smallness 
of the intergenerational mixing elements (5's) when compared with the diagonal elements. 
Otherwise, second or even higher order mass insertions have to be taken into account, as is 
the case for Kaon decays j^]. 

The more recent approach, the general mass eigenstate formalism, is an alternative to the 
mass insertion formalism. Its advantage is that it allows for large off-diagonal elements where 
the mass insertion method is no longer applicable. In the general mass eigenstate formalism, 
the mass matrix in Eq. ()3.4|1 (and the similar one in the down-sector) is diagonalized and the 
flavor changing parameters enter into our expressions through the matrix Tql qji defined 
in Eq. ()3.3|) . So, in the rare top decays t cV, the new flavor changing neutral currents 
show themselves in both gluino-squark-quark and neutralino-squark-quark couplings in the 
up-type squark loops and in the chargino-squark-quark coupling in the down-type squark 
loop. 

After giving expressions for the effective one-loop tcV vertex in the next section, we anal- 
yse the flavor- diagonal case first. Then we consider the effect of intergenerational mixings 
on the rate of the process t — > cV. We assume significant mixing between the second and 
third generations only, in both the up and down squark mass matrices. 
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(a) . (b) 



t 

9,7, z 



(d) 

"k(dk) 



* 1 ► 1 » ► 1 » —7 > 1 



^■*'uh(dh) ^^UkCdk) 

t * ' ' c t 

9,7, Z 

FIG. 1: The one-loop Feynman Diagrams contributing to t — > cg,^,Z including gluino, chargino 
and neutralino loops. For each decay mode, not all four types of diagrams contribute. In t ^ eg 
case, diagram (d) for the chargino and neutralino loops doesn't contribute. In t — > 07 case, diagram 
(d) for the gluino and neutralino loops doesn't contribute. In t — > cZ case, diagram (d) for the 
gluino loop doesn't contribute. 

IV. t ^ IN LRSUSY 

In this section we present the one-loop LRSUSY effective tcV vertex by considering gluino, 
chargino and neutralino loops with up- and down- squarks, including the effects from left- 
right and intergenerational squark mixing. The relevant one-loop Feynman diagrams for 
t cV {V = g, 7, Z) are shown in Fig. ^ Since there are no flavor violating tcV couplings 
in the Lagrangian, there is no tree level contribution. Under the assumption that the charm 
quark mass is negligible, the transition amplitude of the decay is given by 

M{t^ cV) = u{p - k)A'^u{p)e^{k, A), (4.1) 

where p and k are the momenta of the top-quark and gauge boson V, respectively; e{k, A) 
is the polarization vector of the gauge boson. The effective vertex Ay can be expressed as 

K = ^ E h''iPLF^L + PRFvR) + — iPLTl^L + PRTvR)), (4-2) 

where Pl,r = (1 =F 75) /2 is the left (right) chirality projection operator, Fy ^.(r) and Ty ^.(i?) 
are the form factors. In order to compute these form factors, we need to diagonalize the mass 
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matrices for both chargino and neutralino sectors, which are given in Appendix ^ Certain 
Feynman rules of the relevant vertices are presented in Appendix El The calculations are 
carried out in the 't Hooft-Feynman gauge in d-dimensions. Below we give the explicit 
expressions of the form factors for each decay. We prefer to use the notation in Ref. ^ 



which makes the cancellation of infinities more obvious when we set (1 = 4. 
factors into compact form we have further used the identities given in Ref. 

A. t ^ eg decay 

1. Gluino contribution 



?o put the form 

0. 



3^ 87r2 



k=l 



m?, 



p9 



k=\ ^ 



rpg 



C2[G] 
2 

Stt" 



■pfc3^-p*fc^ I" _(jcg[kk] _|_ (jdg[k] _|_ ^[fc] _|_ Qg[k] _|_ ^glk] 



b r 

En \^?^('rk'^ p*fc2^[fcfc] -pkz -p*k2^g[kk] 



k=l 

2 I ^ UL^ UL \ ^gt ~ ^gtov ~ ^g ~ ) ^ ^ UR^ UL\^t + "-^i 



rpg 

^gR 



,g[kk] 
top 



k=l 



n FpWf^S p*fc2^[fcfc] pA;3 p*fe2^9[fcfc] 



+ 



^ UR^ UR ^qt ~ ^ntov ~ ^(1 ~ + ^ UL^ UR\^t + '-^i 



m-, 



•'top 



(4.3) 



2. Chargino contribution 



^gL 



9sg 



167r2 



EE 

k=l a=l 



k,2] 



H 



,,fc,2l 



DR 



Q[a,k,3] _ jy[a,k,3]^ ^ ^cx^[akk] 



'DR 



'-'BE 



[afc]' 
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G 



a,fc,3] 
DR 



H 



a,fe,3] 
DL 



'-'SEG 



'-^SEG V'h 



^gR 



„ „2rpa/ 6 5 
gsQ -t \ - \ ^ f^*[a,k,2] 
- IQ^2 2^2^[^DR 



H 



*[a,k,2] 
DL 



k=l a=l 



I ^[a,k,3] Tj[a,k,3] 
~ I ^DL ~ ^DR 

Omni 6 5 

gs9 T"' 



0) 

[a,k,3] 



mr 



rrit 



Qia,K,3i _ j^[a,k,3]\ f^cx+lakk] _ ^X+K]\ 
DR DL I \ e SE I 



SEG '-^Qpn [lit 



SEG \''''t 



0) 



m" 



mt 



f ^*[a,k,2] rT*[a,k,2]\ ( ^[a,k,3] o-[a,fc,3]\ ^[a. 
167r2 l^l^\yDR -^DL ) [<~-DL '^DR j <^x^ 
k=l a=l L 



kk] 
X+top 



+ [^DR - ^DL ) 



rr[a,fc,3]\ ^x^l^kk] 
^DL j '^top 



rpx+ _ 9s9 J- \- \- ^*[a,k,2] 
k=l a=l 



TT*[a,k,2\\ 
^DR J 



^[a,k,3] 
^DR 



^DL ) 



.kk] 
X+top 



+ ( Ggf'^l 



TT[a,k,3]\ ^x+Kfe] 
^DR J '-^top 



(4.4) 



3. Neutralino contribution 



^x° _ _9s9'T- ^ ^ 

^gL- " 



-i[n,k,3] , TT[n,k,3]\ 
^UL "T -"f/ii ) 



2rpal 6 9 r 

fc=l n=l L 

2rpal 6 9 r 

fe=l n=l L 

^X° _ ( ,/^n*[n,k,2\ TT*[n,k,2] 

9L- 2^2^\y^^UR -^UL 

k=l n=l 

k,3] 



m 



X" 



mt 



A /^/--([",fe,3] I TT[n,k,3] \ ^[nkk] 
\y ^'^UL + ^UR J *-"x° 



X°top 



n -1 



r7-[n,fe,3]\ ^x°["fefe] 
-"[/L J ^top 



2rpal 6 9 p 

= ^ E E (v^^?*J"/'^' + ^S'^^'^O {^^^'^-'' 

k=l n=l L 



_ jj[n,k,3]\ ^[nkk] 
UL j x^top 



(4.5) 
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B. t —>■ cj decay 



1. Gluino contribution 



glgQuC2[F]smd- 



w 



87r2 



k=l 



pfc3 -p*k2 
J- UL^ UR 



, T^k3 p*k2 ( _ fi'gWi /_2 _ n.\\ 
^ ^ UR'- UR I *^SEG ^SEGv'n —^) J 



mt 



mt 
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rpg _ 9l9QuC2[F] sin^w /p^a p*fc2^9[*;fc] pfcs p*fe2^[fcfc] \ 

-'tL — \ f^-f^ UL^top ^ UL^ UL^gtop ) 

k=l 

rpg _ 9s9QuC2[F] sm 6'w /p^a p*fe2^5[fefc] pits p*fc2^[fcfc] 

-'-^R— 2-^ \^ UL'- UR^top UR*- UR^gtoi 



2. Chargino contribution 



■^jR 



sin 



6 5 



k=l a=l 



a,k,2] 
DR 



X 



(5 (7g^^"^' + [aafc] _|_ (^[aafe] _j_ (jX'^lao.k] _|_ ^x+[aa*:]^ _|_ ^Q[o.,k,3] _ ^[a,fc,3]^ 

c^-f 1^ - Qu (c&r - c-s^Er = 0)) ^ 



g'^ sin 



6 5 



167r2 



fe=l a=l 



.kk] 



—Q Cg^'"*' + [""^1 -)- (^H'^l _|_ (^X+[aafcl _|_ (^X+[aafc]^ _|_ ^^Kfc,3] _ ^[a,fc,3]^ 



X ( ^^^h"';^ - Q« {c£t' - C£t\m^t = 0)) 



5''^ sin 6' 



6 5 



167r2 



k=l a=l 



kk] 
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5. Neutralino contribution 
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2. Chargino contribution 



^ ZL 



ZR 



9" 



6 5 



IGtt^ cos^w 



fe=l a=l 



Y*[a,fc,2] 7_7-*[a,fe,2] 
'DL ~ ^DR 



G 



a,A;,3] 



DR 



SE 



h=l 

5 



6=1 ^ 

I fri\P,k,3] TT[b,k,3]\ I /^'[aak]^ab , r^'[bbk]^ab , /-<[ab*;]/n 



a6 



i b 



^7^ 



6 5 



IGtt^ cos 6'w 



jfc=l a=l 



EE(^?/''-^: 



:[a,fc,2] 



G 



a,k,3] 
DR 



_ TT[a,k,3]\ 



0( 



X+[ak] 
SE 



+ (^<^L»L - ^DR ) y-'SEG - ^SEG l"^t - ^) ) J V« Sill ^w 

+ E (^&« " - ^ot'^) ^e"""^""'^ E ^*d1^dl - QdS'' sin^ ^w 

h=l \ j=l y 

5 



6=1 ^ 



15 



^ ZR 



^DL ^DR 



) 



9" 



r 6 



167r^ cos 



/ J / J \ DR DL J / J \ \ DL DR J -^^^-f 



hk] 
X+top 



fc=l a=l 



■ h=l 



6=1 



'DR 



[6,fe,3] 



H 



[6,fc,3] 
DR 



/ ^/[aafc] ^[aafc] ^[abfc] '^t \ 



x"^ 



r 6 



IGtt^ cos 



/ Q*[a,k,2] _ ^*[a,k,2\\ ST^ f ( Qi^XA _ ^[a,h,d.\ \ ^[w 
/ J / J \ DL DR J / J \ \ DR DL j x"* 



hk] 
X+top 



k=l a=l 



■ h=l 



6=1 



'DL ^DR 



^DR ^DL 



^i[bbk]^ab 



I /^l[cLa,k\ ^[aak\ ^[abk\ 

+ I <^^+t + <^^+t - ^x+ 



o 



ah 



(4.10) 



3. Neutralino contribution 
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Here e = 2 — dji id is the number of dimensions), C2[-F] = 4/3(C2[G'] = 3) is the 
quadratic Casimir operator of the fundamental (adjoint) representation of S'[/(3)c with 
iT°-'T°-')ij = C2[F]6ij (T"' are the SU{2) generators in the fundamental representa- 
tion with normalization Tr ^T°-'T^'^ = 5"""/2). The indices a,b run over the intermediate 
chargino mass eigenstates from 1 to 5 while n, m run over the intermediate neutralino mass 
eigenstates from 1 to 9. The indices fc, h are used to represent squarks mass eigenstates, 
which run from 1 to 6. The arbitrary parameter A is introduced in Eq. ()4.12|1 to make the 
argument of the logarithm function dimensionless and our results are independent of A. For 
w = g, throughout the Eqs. ()4.3p - ()4.13|) . we dropped the indices a, b or n,m carried by the 
C functions defined above with a generic w dependency, as gluinos are assumed to have 
identical masses. 

There are two constraints on the form factors that one can check. The first one is the 
cancellations of infinities when we set c? — > 4 (e ^ 0) limit and the second one is the gauge 
invar iance requirement that the coefficient of term in the t ^ cg^cy decays should vanish 
in the /c^ — limit. Obviously this is not the case for the t ^ cZ mode. 
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Cancellation of divergencies can be realized in two ways in d-dimension. Either one can 
adopt the on-shell renormalization scheme (see for example where there is no contribu- 
tion from self energy diagrams (the diagrams (a) and (b) in Fig. and find counterterms 
which exactly cancel the 1/e terms, or include all possible diagrams contributing to the 
process and expect to have overall cancellation of 1/e. We followed the latter. The explicit 
e dependencies are included inside the C-functions defined in Eq. ()4.12j) . As mentioned 
before, this way of writing the form factors has an advantage over other formulations, such 
as Passarino-Veltman, since it makes the cancellation of 1/e terms more apparent. In the 
expressions given above, there are, for example, some combinations of C-functions appear- 
ing repeatedly in which cancellation of 1/e terms requires no computation at all. Note that 
by using the Gordon identity one can decompose Eq. (j4.2|) into an alternative form which 
contains the Dirac structures 7^ and a'^'^k^, where cr^'^ = (i/2)[7^, '-/"]. So, there will be a 7^^ 
component coming from 2p^ part and while we are checking the cancellation of divergencies, 
the combination Fy^j^^ + Ty ^^^^-^ should be considered. However, none of the T'y ^(j^) form 
factors has 1/e dependency so that this point is going to be important only in the discus- 
sion of gauge invariance. Cancellation of 1/e happens in a nontrivial way but is relatively 
straightforward to check for the t eg and t ^ cy decay modes with respect to the case 
for the t — > cZ decay, which requires use of some further properties of Tql,qr, 

3 

E/-p*hi-pki 1 -p*hi-pki \ r/i 
QL^ QL + J- QR^ QR) — " 



i=l 



T.^*qIqR^QL,qr = S'', (4.14) 

k=l 

with unitarity properties of the {U*,V) and A^' matrices defined in Appendix 1X1 in the 
chargino and neutralino sectors, respectively. 

Gauge invariance requires vanishing the 7^^ terms for massless gauge bosons. Thus, we 
should expect to have Fy ^^-^ + T^ri^l) ~ '^'^Y ^ = S'^T in each w = g,x~^yX'^ cases. 
Using the identities given in Ref. jSj for w = g together with the similar ones for w = 
and x° the vanishing of Fy^^j^^ + Ty ^f^^^ for each w is straightforward. 

The branching ratio of the decay t ^ cV can be given, neglecting the mass of the charm 
quark 

\2 
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where = m%/mf and rt()f: — bW) is taken as the total decay width of the top quark. 



(4.15) 



V. NUMERICAL RESULTS 



In this section we present some numerical results for the BR of the decays t — > 
cV {V = g,'y, Z) including the contributions from gluino-up-squark, chargino-down squark 
and neutralino-up-squark loops within the context of both flavor diagonal and non-diagonal 
scenarios. We prefer to discuss each contributions separately, as this is missing from 
the literature. Within the framework we have chosen, there is a large set of parame- 
ters that need to be fixed. Throughout our numerical computations, the SM parame- 
ters are taken as ruf, = 4.5 GeV, = 173.5 GeV, = 80.425 GeV, mz = 91.187 
GeV, as{mz) = 0.1172 and sin^ = 0.2312. The elements of CKM matrix are taken 
as irg^KM = (0.973, 0.220, 3.67 X 10-3) , ir^j^j^ = (-0.224, 0.986, 41.3 x lO^^) , i^t^j^j^ = 
(0.01,-0.05,0.997). The free SUSY parameters in the LR symmetric framework^ are 
Mr, gn, VAji and in addition to the usual ones. Ml, My, Mgusy, A, /i, and tan/5 (we as- 
sume bilinear and trilinear scalar couplings in the soft symmetry breaking Lagrangian to be 
flavor diagonal, fiij = fi, = A). The mass of the gluino needs to be fixed as well. In addi- 
tion to these, there are new flavor changing parameters, (5^[^))jj, iy^j = 2, 3, A, B = L, R 
appearing in both the up-type and down-type squark mass squared matrices defined in 
Eq. ()3.4j) . These are relevant only in the flavor non-diagonal case and we vary them in (0, 1) 
interval. 

The 5x5 chargino and 9x9 neutralino mass matrices are diagonalized numerically and 
the following relations are assumed = Mr, My = ML/2,gji = gi and we set M^ = 150 
GeV, = vsji = 1 TeV, /i = 200 GeV, and tan/3 = 10 throughout our analysis.^ There 

The LRSUSY is assumed not embedded into some larger supersymmetric grand unified theories hke 
50(10), Eq. Otherwise, one can relate some of these parameters by using the relations among them at 
the unification scale with the use of renormalization group equations. Such framework, however, would 

lead to more complicated particle spectra. 
^ We also considered large tan/3 values and discuss it at the end of Section IV Bl 
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are experimental lower bounds on the masses of charginos, neutralinos as well as gluinos and 
squarks j^]. For the values of the parameters given above, the lightest chargino is around 
130 GeV and the lightest neutralino 90 GeV, which are consistent with the experimental 
lower bounds 2^. To determine the squark masses, we need to further fix Mgusy and the 
coupling A. We set A = Mgusy and vary them within (100, 1000) GeV in the fiavor diagonal 
case, and set Mgusy = 300,400,1000 GeV in the fiavor non-diagonal case. For example, 
the lightest up-type (down-type) squark mass is around 250 (290) GeV for Mgusy = 300 
GeV in the fiavor diagonal context. When the fiavor changing parameters are turned on, 
the lightest up- type (down-type) squark is around 100(200) GeV for intermediate values of 
^RR,LL _ ^RL,LR^^RR,LL _ ^RL,LRy y^y^qyi the fiavor violation comes only from the RR or 

LL sector, we are getting slightly larger squark mass values. We now discuss in the following 
subsections the fiavor diagonal and non-diagonal cases separately. 



A. Flavor Diagonal Case 

We discuss first the fiavor diagonal case and analyse the fiavor changing effects based 
on the fiavor diagonal results. We assume that the SM contribution to the t cV {V = 
g, 7, Z) decays is negligible with respect to the contributions from the LRSUSY as the 
total branching ratio is concerned. In some cases, chargino or neutralino loop contributions 
might be comparable with the SM values. However, we mainly concentrate on the relative 
contributions from each loop involving supersymmetric particles. In all our considerations, 
the gluino loop dominates both chargino and neutralino loops. 

We first investigate the t ^ eg branching ratio as a function of Mgusy within the range 
(100, 1000) GeV for tan/3 = 10 and nig = 300 GeV. Fig. |2 shows the dependency including 
the gluino, chargino and neutralino contributions separately. Neutralino dominates the 
chargino contribution for the most part of the parameter space (except for large Mgusy values, 
where they are roughly equal). There is a two orders of magnitude difference between them 
for Msusy ~ 200 GeV which gets smaller as Mg^sy gets larger. The gluino contribution is 
around three orders of magnitude larger than the neutralino contribution for Msusy ~ 200 
and the difference gets even bigger, reaching 4 x 10~^ when Msusy is around 100 GeV. 
Overall both the neutralino and chargino contributions are practically negligible. In general, 
the neutralino contribution is bigger than the chargino contribution in the flavor diagonal 
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FIG. 2: The gluino, chargino and neutrahno contributions to the BR(t — > eg) as a function of 
Msusy in the flavor diagonal scenario for rUg = 300 GeV, tan/3 = 10 together with the parameter 
values defined at the beginning of section The solid, dashed, dotted-dashed curves as coded 
represent gluino, neutralino and chargino contributions, respectively. 

scenario for t —>■ cV (with the notable exception of the t cZ decay mode and the t —>■ cy 
decay for certain Mgusy values) and this is consistent with the findings of Lopez et al. P|. 
Note that this picture is reversed for the rare decays involving down-type quarks like bottom- 
strange quark transitions (6 decays). There the chargino contribution is always larger than 
the corresponding neutralino one. In h decays, the chargino couples with up-type squarks 
while here it couples with down-type squarks, whereas the neutralino does the opposite. So, 
;he features seen here are consistent with the characteristics of the down-type rare decays 

In Fig. El we present the Mgusy dependency of the branching ratio of the decay t — > 07 
for the same parameter values chosen above. The gluino contribution is one to two orders of 
magnitude smaller than the one in the t eg decay and it reaches 2 x 10"^ at the maximum 
level. The neutralino contribution is two to four orders of magnitude smaller than gluino. 
Even though, after Msusy ~ 300 GeV, the chargino contribution is larger than the neutralino, 
both are still negligible in the entire parameter space considered. 

The final figure of this subsection. Fig. HI shows the same dependency, but for t cZ de- 
cay. The general pattern is the same except that chargino dominates neutralino everywhere 
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FIG. 3: The same as Fig. IJJbut for the t — > 07 decay mode. 

in the (100, 1000) GeV range for Mgusy (one order of magnitude larger than the neutrahno) 
and furthermore the chargino loop contribution is almost constant as Mgusy changes. The 
gluino contribution takes maximum values of 7 x 10"^ but decreases sharply as Mgusy gets 
bigger and even becomes comparable to the chargino contribution. 

When we consider the total branching ratios of the three decay modes, in the flavor 
diagonal case, the largest one is t ^ eg as expected and the smallest is the t ^ cj decay 
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FIG. 4: The same as Fig. |21but for the t ^ decay mode. 
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mode. For instance, we have BR{t eg) ~ 10"^ BR(t cZ) ~ 5 x 10-^ BR{t 07) ~ 
10^^ for nig = 300 GeV and the intermediate Msusy = 300 GeV value . It is also possible to 
get one to two orders of magnitude larger values for some smaller gluino mass. In all three 
cases, both the chargino and neutralino contributions are quite suppressed and negligible 
for practical purposes. 

Before discussing the flavor non-diagonal effects, we would like to comment on the current 
experimental limits on these decay modes. The best bound for t cy decay mode is from 
CDF which looked for flavor changing top quark interactions in pp at 1.8 TeV center 
of mass energy. The bound is BR(t 07) < 3.2 x 10^^ (95%CL). The bound for the 
t cZ decay channel is weaker, BR(t cZ) < 0.137 (95%CL) for rrit = 174 GeV, from 
the OPAL experiment [2^ which searches for single top quark production in the e^e~ —>■ tc 
reaction at around 200 GeV center of mass energy.^ For the top quark decaying into the 
charm quark and gluon, there is no current experimental bound available since background 
problems make detection of such channel difficult l30l l. 

These bounds will hopefully be improved by LHC in near future m^. Both LHC and the 
future LC have an advantage over some other colliders for detecting rare top decays because 
of having better statistics capabilities with lower backgrounds. For example, the foreseen 
sensitivities to these channels can be as small as 30|, |31 1 

BR{t eg) < 7.4 x lO^^ , 
BR(t C7) < 1.0 X 10"^ , 

BR(t cZ) < 1.1 X 10"\ (5.1) 

so that consideration of the flavor non-diagonal effects are essential, since the enhancement 
from such effects makes it possible to reach the experimentally accessible range for at least 
some of the modes considered here. Theoretical estimates for top flavor changing neutral 

n 

couplings have been performed by using parton- level simulations [24. Further references 
and a description of the processes involved can be found in 3^. 



There is a slightly weaker bound from the ALEPH experiment 
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B. Flavor Non-Universal Case 

he remaining part of our study is devoted to discussing the effects of flavor changing 
mixings in both the up- and down-type squark mass matrices on the rare t ^ cV {V = 
g, 7, Z) decays. As motivated in section IllH we only concentrate on mixings between the 
second and third generations and neglect any kind of mixing involving the first generation. 
Furthermore, unlike some previous studies where the so-called "mass insertion" method has 
been used, we follow the general "mass eigenstate" formalism. In addition to the parameters 
in the flavor diagonal case, we have eight more, essentially unknown, parameters in each 
sector as given in Eq. (j3.4j) (that is, there are totally sixteen parameters in both the up and 
down sectors). 

In order to reduce the set, we should conservatively make some further assumptions and 
consider some limiting cases. First of all, we set the flavor and chirality diagonal elements 
of the squark mass matrix to a common SUSY scale, Mgusy, which enables us to define 
dimensionless flavor changing parameters, {5Q^)ij, i^j = 2,3, A, B = L,R. In addition 
to that, we are going to restrict our consideration to a LRSUSY model with Hermitian 
or symmetric trilinear couplings so that the previously defined mixing parameters satisfy 
i^Q^)ij ~ i^Q^)jiy "^T^i = 2,3, A,B = L,R. Furthermore, from the LR symmetry we 
expect {SQ^)ij = {SQ^^)ij to hold, though we will not always assume it. Note that this is in 
general not necessarily true (MSSM is an example in this respect). Therefore, we are left 
with four additional parameters, (5^^)23, (5^'^)23, {^1/^)23, {^1/^)23, in the up squark sector, 
and similarly for the down squark sector. 

We investigate five limiting cases: the case where the dominant mixing effects come from 
only the LL, or RR, or LL+LR, or RR+RL, or LR+RL blocks of the squark matrices. 
However, only the RR, RR+RL and LR+RL cases are presented here^ since the ones with 
LL are very similar but not identical. This is basically because the flavor conserving LR 
mixing in the second generation is not equal to the one in the third generation. We don't 
consider limiting cases between 5's and flavor conserving but chirality changing entries of the 
matrix in Eq. (j3.4j) . One reason is that we kept the trilinear soft terms set to A = Msusy which 
makes such limiting cases very similar to the flavor diagonal case which we have discussed 

^ We concentrate on mainly the RR sector since the model has an additional SU{2)ii group, which makes 
the right-handed sector more interesting. 
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in the previous section. So, we proceed to investigate the five hmiting cases in each decay. 
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FIG. 5: The gluino, chargino and neutralino contributions to the BR(t eg) as function of 
('^l^(D))23 case where the flavor changing effects come from the mixing in the RR sector 

only. The graph on the left(right) denotes the gluino(neutralino) contribution for the up RR 
mixing. The one in the middle is for the chargino as a function of (5|)^)23- Three representative 
Msusy values (300,400,1000) GeV are chosen with m-g = 300 GeV, tan/3 = 10, vl = M^usy and 
/i = 200 GeV. 



In Fig.0 we plot in a set of three graphs the gluino, chargino and neutralino contributions 
to the BR(t eg) as a function of (5/7^)23, for Msusy = 300,400,1000 GeV, assuming all 
other mixings to be zero; nig = 300 GeV and tan/? = 10 are assumed. As seen from 
the graph on the left, the gluino contribution to the branching ratio depends strongly on 
{^u^)23 and gets enhanced as (5^^)23 takes larger values. When we compare with the fiavor 
diagonal case (Fig.|21), for Msusy = 300 GeV, there is at least a factor of five enhancement for 
intermediate values of {6^^)23 and of one order of magnitude for larger values. As Msusy gets 
larger, the enhancement is two to three orders of magnitude for intermediate-large values of 
(5^'^) 23. Two comments are in order. The first one is that the comparison with the fiavor 
diagonal case is not completely right. Here we assume all diagonal elements in the main 
diagonal identical (to Msusy), which leads to larger values for the branching ratio. That is, 
the enhancement would be slightly different if we relaxed this condition. The second one is 
that the full (0,1) interval chosen for the fiavor changing parameters is not always available 
since squark masses become unphysical when 5's exceed certain critical values, which also 
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depend on Mgusy. 

In the same figure, for tlie same parameter values, we depict tlie cfiargino and neutralino 
contributions to tlie BR(t eg) in tlie second and third graphs, respectively. The chargino 
contribution is not very sensitive to 6 for {6^^)23 < 0.8 and the enhancement with respect to 
the flavor diagonal case is one order of magnitude in that range but becomes three orders of 
magnitude after. When we compare it with the gluino contribution, it remains suppressed 
even under extreme conditions (i.e., {S§^)23 = 0, {6^^)23 ~ 1). So, the flavor mixing effects 
from down squark sector do not contribute to the total BR of the t ^ eg decay nearly as 
much as the ones from the up squark sector. For the neutralino case, even though there is 
three to four orders of magnitude enhancement for intermediate or large (^1^^)23 values, it is 
at least one order of magnitude smaller than the gluino contribution. Neutralino contribution 
still dominates the chargino one unless (5^^)23 = 0, {6^^)23 ~ 1 occurs. 
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FIG. 6: The same as FigElbut for the case where both (^^(^^23 and {^y\p^^23 contribute, with 
the assumption (<5^p))23 = (5^(b))23- 



In Fig.ini we show the same as Fig.^lbut this time additionally the mixing in the LR sector 
is turned on. For simplicity we present the case (5^(|)))23 = (5^(£())23- With respect to the 
case where only the RR mixing is turned on, except the chargino contribution, the gluino 
and neutralino contributions are suppressed about one order of magnitude. However, in 
the down sector, the chargino contribution depends strongly on (5|5"^)23 and reaches 10~^ at 
maximum level. In the rest of the interval, it is still very suppressed in the total BR(i(: eg). 



The last case for the BR(t eg), as a function of (5^(^))23 = {SijfD))'^^' depicted 
in Fig. 13 We have smaller contributions with respect to both the RR and RR+RL cases 
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FIG. 7: The same as FigElbut for the case where both (5^|^-))23 and {6^^jj^)23 contribute, with 
the assumption ((5^(|)))23 = (5^fz)))23- 

except for the neutralino contribution, which reaches values obtained in the RR+RL case 
for upper values of {6^^)23 at different Msusy values. So, overall, the gluino gives the largest 
contribution among three in all limiting cases considered above, and the total BR can get 
up to a few times 10~^ for rrig = 300 GeV. 



2. t^c-f 
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FIG. 8: The gluino, chargino and neutralino contributions to the BR(i cj) as function of 
((^^^^23 for the case where the flavor changing effects come from the mixing in the RR sector 
only. The graph on the left(right) denotes the gluino(neutralino) contribution for the up RR 
mixing. The one in the middle is for the chargino as a function of {6^^)23- Three representative 
Msusy values (300,400,1000) GeV are chosen, with rUg = 300 GeV, tan/3 = 10, ^ = Msusy and 
/i = 200 GeV. 

We analyse the t ^ cj decay in the same order as we did the t ^ eg decay and we present 
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only the RR and RR+RL cases. In Fig. |S1 we show the dependence of the gluino, chargino 
and neutrahno contributions in the BR(t cy) to the parameter (5^(^-))23 at different 
Msusy values, 300,400, 1000 GeV for tan/? = 10. The gluino contribution is now one order 
of magnitude bigger than in the flavor diagonal case for Msusy = 300 GeV and {6^^)23 ~ 0.5 
and it can go up to three orders larger for Mgusy = 1000 GeV with maximum possible mixing 
in the up RR sector, or for smaller Mgusy with slightly smaller mixings, when it reaches 10~^. 
In the chargino case, in the down sector, the dependence of the BR to Mgusy is very weak 
and enhancement with respect to the flavor-diagonal case varies from two to three times to 
about two orders of magnitude. Here the curve for Msusy = 300 GeV has a minimum which 
is peculiar only to this Msusy value and occurs due to some precise cancellations in the loop 
functions. For example, we don't have such behavior in the down LL mixing case which 
we haven't discussed here explicitly. It is, however, possible to get bigger enhancements for 
smaller Msusy values. For the neutralino contribution, bigger enhancements ranging between 
two to four orders of magnitude occur, depending on the value of Msusy and the amount 
of mixing allowed in the up RR sector. Unlike the flavor diagonal case, the neutralino 
contribution always dominates the chargino one through the entire interval scanned here; 
while it is still suppressed with respect to the gluino contribution. 




FIG. 9: The same as FigElbut for the case where both (5^^^)23 and (5^^))23 contribute, with 
the assumption ((5^^p23 = (5^(b))23- 

The RR+RL case (((5^(^))23 = ('^t/(i)))23) ^^e BR(t 07) is similar and shown in 
Fig. 121 The gluino contribution remains always one order smaller than in the RR mixing 
case. This is true for neutralino as well even though there is a sharp dependency on the 
mixing parameters in the up sector. In this case, like in the flavor diagonal case, the chargino 
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curves can still cross the neutralino curves and become larger with larger down type RR+RL 
mixing. If one considers the mixing in the up and the one in the down sector to be completely 
independent, one could end up with different conclusions. 



3. t^cZ 
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FIG. 10: The gluino, chargino and neutralino contributions to the BR(t cZ) as function of 
('^C/(D))23 case where the flavor changing effects come from the mixing in the RR sector 

only. The graph on the left(right) denotes the gluino(neutralino) contribution for the up RR 
mixing. The one in the middle is for the chargino as a function of {8^^)2'i- Three representative 



M, 



susy 



values (300,400,1000) GeV are chosen with m-g = 300 GeV, tan/3 = 10, ^ = M^usy and 



^ = 200 GeV. 

Finally, we present the final decay channel, t cZ. Fig. ^1 shows the branching ratio 
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FIG. 11: The same as Fig^Jbut for the case where both ((5^^^)23 and ((5^^p23 contribute, with 
the assumption {5^fj~,-^)2z = {^uiD))'^^- 
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of t — s> cZ as a function of ('5^(|)))23 when the rest of flavor changing parameters are set 
to be zero. The gluino contribution becomes enhanced one to more than four orders of 
magnitude compared to the case where all flavor changing parameters are turned off. A 
10~^ branching ratio seems to be reachable. We can the discuss neutralino contributions 
and directly compare with the gluino contributions since, unlike the chargino contribution, 
they vary with respect to the same mixing parameter. The neutralino contribution is more 
than two orders of magnitude enhanced over the flavor diagonal case but still suppressed 
with respect to the gluino contribution. The effect in the chargino sector is tiny. 

The RR+RL case, presented in Fig. ^2 is very similar to the above case for both the 
gluino and neutralino contributions but the mixing in the RR+RL down sector contributes 
to the chargino loop more effectively and increases its contribution to one order of magnitude 
with respect to the RR mixing case considered above. The chargino contribution, like in the 
flavor diagonal case, is still dominant over the neutralino one. 

Consideration of flavor changing effects in either sectors is quite important for rare top 
decays t —>■ cV {V = g,'~f, Z) not only for the expected enhancements, but also for revealing 
the details of the SUSY flavor changing mechanism. 

Before closing this section, we would like to comment on the large tan/3 effects in both 
flavor diagonal and non-diagonal scenarios. As can be seen from the Feynman rules for 
gluino, the gluino loop for each decay is not sensitive to tan (3 value. Basically, the branch- 
ing ratios are slightly bigger for tan/3 values smaller than 10 and almost constant after 10. 
This indeed makes the total branching ratios of the decays insensitive to tan/5. Since we 
discuss contributions from each loop separately, it is worth mentioning the tan (3 dependen- 
cies of chargino and neutralino loops. We comment on specifically the t ^ cy decay, but 
these comments apply to other decay modes as well. In the flavor diagonal scenario, unlike 
neutralino case, tan/5 dependency of the chargino contribution to t — > 07 is quite strong. 
If one scans the BR(t C7) in tan/5 in the — 50 interval, the chargino contribution can 
reach 2 x 10-^^{3 x 10"") for Msusy = 300, 400 (1000) GeV, which is about 40 (30) times 
bigger than tan /5 = 10 case and three orders of magnitude bigger with respect to very small 
tan (3 values. For the neutralino case, the branching ratio of t — > 07 is almost constant as 
tan/5 changes for Mgusy = 300,400 GeV. As Mgusy gets bigger, there is a slight increase in 
the branching ratio as tan j3 increases. So, for Msusy = 1000 GeV, we have one order of mag- 
nitude enhancement at maximum value of the branching ratio in the interval tan j3 G (0, 50) 
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and five times bigger branching ratios when we compare the case for tan (3 — 50 with the one 
for tan/3 = 10. Even though the chargino contribution becomes larger for large tan/3, it is 
still quite suppressed with respect to the gluino contribution. For example, for tan /3 = 50, 
the chargino contribution is around 75; 40; 4 times smaller for Mgusy = 300; 400; 1000 GeV. 
Where the chargino contribution was one to two orders of magnitude smaller than the neu- 
tralino contribution for small tan/3 values (tan/3 < 5), it becomes more than one (two) 
order(s) bigger then the neutralino contribution for Mgusy = 300, 400 (1000) GeV as tan/3 
takes its maximum value considered here. 

For the fiavor non-diagonal case, we consider again as an example the t ^ decay 
with mixing only in the RR sector. The gluino contribution is very similar to the flavor 
diagonal case. The chargino contribution is again most sensitive to tan/3 among the three 
contributions. The sensitivity in this case is slightly less, because we have assumed for 
simplicity the diagonal entries of the squark mass matrices identical to Mgusy For the fiavor 
mixing parameter {6^^)23 — 0.5, the chargino contribution to the branching ratio becomes 
35; 25; 10 times bigger for Mgusy = 300; 400; 1000 GeV when we compare the cases tan P — 10 
and tan /3 = 50. For the neutralino case, the discussion is very similar to the fiavor diagonal 
case with a reduced dependency on tan/3. In summary, when we add up all contributions 
the branching ratio is not sensitive to the tan/3 parameter because the gluino loop which 
dominates all contributions and it is not sensitive to tan /3. When we look at the individual 
contributions to the branching ratio, the chargino contribution has the strongest dependency 
on tan/3 both in fiavor diagonal and non-diagonal scenarios. Thus, our overall analysis is 
valid as one varies tan /3 in the phenomenologically relevant interval. 

VI. CONCLUSION 

We have presented a complete analysis of the two-body fiavor violating decays of the top 
quark in a fully left-right supersymmetric model. The model has a left-right scalar quark 
sector, as well as a right-handed gaugino in both neutral and charged sectors. 

We have first evaluated the branching ratios in the fiavor diagonal case, which corresponds 

to the constrained LRSUSY, where explicit flavor violating terms do not exist in the squark 
mass matrix, and the only source of flavor violation is the CKM matrix. The enhancement 
due to the right handed sector proves to be minimal, and top quark FCNC decay is somewhat 
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below the detectable level at future collider experiments. We present the results for the 
gluino, chargino and neutralino contributions to the decays t — > eg, 07 and cZ explicitly and 
separately, unlike previous analyses. We confirm that branching ratios can reach 10~^ for 
the gluon, and 10^^ for the photon and Z, in agreement with previous results. 

We then analysed the branching ratios for t — > cV in the unconstrained LRSUSY, where 
flavor changing elements in the squark mass matrix are allowed to be arbitrarily large, 
though for the mixing between the second and third generation only. This mixing is rather 
unconstrained in the top-charm squark sector, and only very weakly constrained in the 
bottom-strange squark sector from b decays. We analyse again separately contributions from 
gluino, chargino and neutralino to t ^ eg, c-f and cZ and look at results for the branching 
ratios when flavor violation is driven by the right-right, right-left and left-right, or both, 
squark mixings. As expected, enhancements of the branching ratio could occur here, and 
we obtain, under the best circumstances, BR{t — > eg) close to 10~^, while BR{t — > 07) and 
BR{t — > eZ) can reach 10~^ and 10~^ respectively. These enhancements are comparable 
to other studies in the unconstrained MSSM, so while t decays are important as a probe 
of SUSY FCNC couplings, they do not appear to differentiate between the MSSM and 
LRSUSY, unless there are clear indications that it is the right-handed scalar quark sector 
which is responsible for flavor violation. In this case LRSUSY provides large enhancements, 
at least competitive with ones driven by the left-handed scalar quark mixings in MSSM. 

Comparing t to b decays, it appears that there the chargino contribution is more evident 
and sometimes comparable to the gluino; while this is never the case for t ^ eV and thus 
here the effect of the right-handed gaugino is relatively obscured. 
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APPENDIX A: CHARGINO AND NEUTRALINO MIXINGS 

In this appendix we present the chargino and neutralino mass matrices and the way of 
their diagonalizations. 
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1. Chargino Mixing 



The terms relevant to the masses of charginos in the Lagrangian are 

+ H.c. , 



1 , I T-i, / X'^ 
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(A2) 



where we have taken, for simphfication, fiij = fi. Note that although there are six charginos, 
and 5j decouple from the spectrum because we have chosen the VEV's of Ax, and 5l to 
be zero. 

The chargino mass eigenstates Xi obtained by 



xt = Vij'^t, Xi = Uijipj , i, J = 1, ... 5, 
with V and U unitary matrices satisfying 

U*XV-^ = Xd, 



(A3) 



(A4) 



The diagonalizing matrices U* and V are obtained by computing the eigenvectors corre- 
sponding to the eigenvalues of XX^ and X^X, respectively, to the eigenvalues of XX^ and 
X^'X, respectively. 

2. Neutralino Mixing 

The terms relevant to the masses of neutralinos in the Lagrangian are 



(A5) 



34 



where ip^'^ = (— , —i^R, —i^ 
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(A6) 

As in the case of charginos, A° and Si decouple from the spectrum because we have chosen 
the VEV's of and 6l to be zero. The mass eigenstates are defined by 



X° = Ar,,V'° (z,j = l,2,...9), 
where is a unitary matrix chosen such that 

NZN^ = Zn, 

and Zd is a diagonal matrix with non-negative entries.^ 

One can further switch to a basis involving the photino, left and right zino states 

= i-zx^, -a^^, K% ~5% 0?„ 
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where the photino and left and right zino states are defined as 
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(A9) 



The positivity of the entries of Zu can be achieved by muhiplying some of the rows of N with i. This 
freedom comes from the fact that in order to determine iV the square of Eq. (|A8I) needs to be considered. 
We have a similar situation for the chargino sector. For example, see 13 for details in the MSSM. 



35 



The mass matrix in the above basis becomes 
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with the functions 
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(All) 

where we have used = and gfy = e/\/ cos 2^w • Furthermore, the unitary matrix iV can 
be expressed in the rotated basis as 

/ N'\ ( sin^w sin^w a/cos 2^w \ / iV,^ \ 
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APPENDIX B: FEYNMAN RULES 



In this appendix, we give the Feynman rules involving SUSY particles. The ones for the 
SM particles can be found in many textbooks. 

• Gauge boson-squark-squark interaction 
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Gauge boson-chargino-chargino interaction 
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• 9iX^} x"'"-quark-squark interaction 
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